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Magnetic susceptibilities ol spin-(S', s) ferrimagnetic Heisenberg chains are numerically investi- 
gated. It is argued how the ferromagnetic and antiferromagnetic features of quantum ferrimagnets 
are exhibited as functions of (S,s). Spin-(5', s) ferrimagnetic chains behave like combinations of 
spin-(S — s) ferromagnetic and spin-(2s) antiferromagnetic chains provided S = 2s. 
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Quantum behavior of mixed-spin systems is one of the 
hot topics and a recent progress in the theoretical un- 
derstanding of it deserves special mention. In contrast 
with pioneering calculations H] in the 1980s, the current 
vigorous argument might more or less be motivated by 
the Haldane conjecture in that it covers not only ferri- 
magnetism but also mixed-spin antiferromagnetism. As 
coexistent spins of different kinds do not necessarily re- 
sult in magnetic ground states, it remains a stimulative 
problem there whether the system is massive or massless. 
Actually various mixed-spin chains |],|| and ladders || 
with singlet ground states have also been discussed ex- 
tensively via the nonlinear-cr-model technique with par- 
ticular emphasis on the competition between massive and 
massless phases. On the other hand, an arbitrary align- 
ment of alternating spins S and s, which is described by 
the Hamiltonian 

H = Y J (-hS ] -s 3 +J 2 s r S ]+1 ), (1) 
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shows ferrimagnetism instead of antiferromagnetism and 
thus has gapless excitations from the ground state. 

Chemical explorations into ferrimagnetic chain com- 
pounds have vigorously been reported so far. In an 
attempt to obtain bimetallic complexes including one- 
dimensional systems, CuMn(S2C 2 02)2(H20) 3 -4.5H 2 
H was first synthesized. Since then numerous bimetal- 
lic chain compounds [Q were systematically synthesized 
and measured. In the series of investigations, the smaller 
spin s was fixed to 1/2 and the larger spin S was changed 
from 1/2 to 5/2 featuring several important concepts. 
One of their main interests was the achievement of a fer- 
romagnetic interaction between nearest-neighbor metal 
ions, namely, the design of molecular-based ferromagnets. 
Therefore quantum ferrimagnetism was not necessarily 
the main subject there. Though there also appeared the- 
oretical investigations |l],|8| in close contact with exper- 
iments, the argument was restricted to the spin com- 
binations (S, 1/2) and the larger spin S was sometimes 
treated as classical. 

In such circumstances, several renewed interests [p|-|l9||, 
have recently been introduced into the field. Alcaraz and 
Malvczzi H predicted the universal low-energy physics 



of anisotropic ferrimagnetic chains and the distinct be- 
havior of isotropic ones. Kuramoto jl6| discussed critical 
properties of the model in a field showing its magnetiza- 
tion curve with a plateau. Maisinger et al. JTs| ] revealed 
the field-induced double-peak structure of the specific 
heat. Another noteworthy subject, which is the main 
interest here, is the coexistence of the ferromagnetic and 
antiferromagnetic features in quantum ferrimagnets. The 
author et al. pointed out in their recent study Jl2|,[l4],[l7| 
on the model (Q) with (S, s) = (1, 1/2) that at low tem- 
peratures the model acts as a ferromagnet, while at in- 
termediate temperatures it behaves like a gapped anti- 
ferromagnet. In this article, we quantitatively verify to 
what extent this scenario is valid discussing various com- 
binations of (5, s). Even though there is an ambiguous 
belief that the (S,s) — (1,1/2) case is representative of 
all the Heisenberg ferrimagnetic chains, it is not yet re- 
ally clear whether or how the quantum behavior of the 
model depends on the constituent spins. We here show 
that the magnetic behavior of spin-(5*, s) quantum ferri- 
magnets is regarded as the combination of those of spin- 
(5 — s) ferromagnets and spin- (2s) antiferromagnets pro- 
vided S — 2s. The dual aspect is qualitatively generic to 
all the Heisenberg ferrimagnets, but those with S = 2s 
should quantitatively be distinguished from the others. 

We investigate the zero-field magnetic susceptibility \ 
of the model ([!]) of N elementary cells with J\ = J 2 = J 
under the periodic boundary condition. We take /jb as 
the Bohr magneton and set the g factors of both spins 
S and s equal to g. We assume that S > s. The nu- 
merical procedure of the quantum Monte Carlo method 
employed has been detailed elsewhere [^l[ . Since the cor- 
relation length of the system is considerably small |l0|,|l2| , 
the thermal quantities in general show no significant size 
dependence. We have calculated the chains of N = 24 
and N — 32, but we find no difference beyond the numer- 
ical uncertainty between them in the temperature region 
treated. On the other hand, for the sake of argument, we 
calculate the susceptibilities of Heisenberg ferromagnetic 
and antiferromagnetic rings as well, where we consider 
chains of N spins. 

We show in Fig. [l] temperature dependences of the 
product xT at various combinations of (S,s). We note 
that here \T shows a minimum in its temperature depen- 
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dence. Although the low-temperature divergence, which 
is proportional to T^ 1 , is reminiscent of the ferromag- 
netic susceptibility, \T approaches its paramagnetic be- 
havior [S(S + 1) + s(s + l)]/3 showing the antiferro- 
magnetic increase. Considering that \T is a monotoni- 
cally decreasing function in ferromagnets, while a mono- 
tonically increasing function in antiferromagnets, the 
present observations may be regarded as a ferromagnetic- 
antiferromagnetic crossover. 

This mixed feature is a consequence of the dual struc- 
ture of excitations |T4| ]. Quantum ferrimagnets exhibit 
two distinct branches of elementary excitations, one of 
which reduces the magnetization M = (S| + s* ) and 
is thus called ferromagnetic, and the other of which en- 
hances that and is thus called antiferromagnetic. The 
lowest-order spin-wave approach Jl(],|l2| to the model (|l|) 
with J\ = J2 = J indeed results in two dispersions 



J 



(S - s) 2 + ASs sin 2 (ak) t(S- s) 



(2) 



with the lattice spacing a, which are quadratic at small 
fc's and gapped from the ground state, respectively. An- 
other naive but suggestive calculation |jlj] is the pertur- 
bation from the decoupled-dimer limit, where the ground 
state and its elementary excitations are trivially obtained 
(Fig. ||). When we turn on the exchange interaction be- 
tween the dimcrs, the localized excitations begin to hop 
to neighboring unit cells and become dispersive. In the 
case of (S, s) — (1,1/2), for instance, setting J\ and J2 
for J and S J, respectively, we obtain 

= - JS[1 - cos(2ak)} + 0(S 2 ) , (3) 
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cj+ = I J + i JS[7 - 12 cos(2a/c)] + 0{5 2 ) . (4) 
2 18 

For an arbitrary combination of (S, s), the similar argu- 
ment ends up with qualitatively the same dispersions. 

Based on this scenario for low-lying excitations, we 
propose in Fig. || an idea of decomposing ferrimagnets 
into ferromagnets and antiferromagnets, where we let the 
decoupled dimers and the Affleck-Kennedy-Lieb-Tasaki 
valence-bond-solid states J20J symbolize ferrimagnets and 
integer-spin gapped antiferromagnets, respectively. Now 
we are led to expect spin- (5, s) ferrimagnets to behave 
like combinations of spin-(2s) antiferromagnets and spin- 
(S — s) ferromagnets. In this context, we compare in 
Fig. ^ the ferrimagnetic susceptibilities with those of fer- 
romagnets and antiferromagnets. Surprisingly, the sum 



of the spin-(2s) antiferromagnetic and spin-(S f — s) ferro- 
magnetic susceptibilities excellently reproduces the spin- 
(S, s) ferrimagnetic behavior on condition that S = 2s, 
whereas such a naive relation does not stand up in the 
other cases. Considering several suggestive calculations, 
we may be convinced of the present observations. 

First, the high-temperature paramagnetic behavior 
distinguishes the S = 2s cases from the others. At high 
enough temperatures, we should obtain the paramagnetic 
susceptibilities as 



fc R TY (S ' s:i ~ fcrri - 1 

NM °sM g + 1 > + '(' + 1 >]' 
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A^ 2 = z(S-s)(S-s+l), 
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(5) 
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It is interesting that we have an equality X s fcrn ' = 

x (S- S )-ferro. + x (2s)-a„tiferro. j ugt j n the cases Q f g = 2s 

Second, the spin-wave dispersions (0) give us a use- 
ful piece of information. The small-momentum ferro- 
magnetic excitations must dominate the low-temperature 
thermodynamics and therefore the curvature of their dis- 
persion is of great interest. Let us recall that spin- S fer- 
romagnetic Heisenberg chains exhibit the dispersion of 
single-magnon excitations, 



ujk = 2JS(\ — cos/ca) . 



(8) 



Defining the curvature v for quadratic dispersions as 
^k^o = v(ak) 2 with the unit-cell length a, we obtain 



(S,s)-feTri. 



Ss 



-J. 
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w (S- S )-fcrro. _ (5_ S )J. 



(9) 
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Here we learn that v (S,s)-fcm. _ w (S-s)-ferro. p rov ided 
S = 2s. Though the coincidence holds within the lowest- 
order spin-wave theory, it looks meaningful. 

Finally, we point out that the modified spin- wave treat- 
ment p2| of spin chains is also suggestive in our argu- 
ment. The modified spin-wave theory was first intro- 
duced by Takahashi [^3| for quantum ferromagnets. Im- 
posing a constraint on the magnetization, he obtained the 
low-temperature susceptibilities of spin-S* ferromagnetic 
chains as 
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+ or i/2 ), 



(ii) 



where t = k^T/J and Q(z) is Riemann's zeta function. Equation ( pd| ) was shown to precisely agree with the ther- 
modynamic Bethe-ansatz calculation 24 in the case of S — 1/2. On the other hand, the author and Fukui p"^ ] 
have recently developed the scheme for quantum ferrimagnets. They constructed the thermodynamics based on the 
dispersions (0) and obtained the expression for the susceptibility, 
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(S,s)-ferri. 
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For ferrimagnets the constraint is not so trivial as for 
ferromagnets. However, as far as the low-temperature 
behavior is concerned, we obtain the universal descrip- 
tion regardless of the artificial constraint. Equation Jl2] ) 
is thus reliable to a certain extent. Now we again en- 
counter the surprising coincidence between the two re- 



sults, x 



(S,s)-ferri. 



x 



(S — s)-ferro. 



provided S = 2s. 



Thus we have learned that several analytic treatments, 
as well as numerical calculations, distinguish ferrimag- 
netic chains with S = 2s from the others. We conclude 
an approximate sum rule among ferrimagnets, ferromag- 
nets, and antiferromagnets on condition that S = 2s: 



X 



(S,s)-fcrri. 



(S— s)-fcrro. . (2s)-antifcrro. 
A i A 



(13) 



Equation (|13j) well holds in a wide region of temperature 
except for the ferromagnetic-to-antiferromagnetic tran- 
sitional period. Both ferromagnetic and antiferromag- 
netic features generically lie in quantum ferrimagnets, to 
be sure, but the quantitative identity of them is estab- 
lished only in the cases of S — 2s. For S > 2s quan- 
tum ferrimagnets are rather ferromagnetic than antifer- 
romagnetic, whereas for S < 2s vice versa. In fact we 
find in Fig. [j] that \T shows a relatively significant low- 
temperature divergence for S > 2s, while it increases rel- 
atively rapidly at intermediate temperatures for S < 2s. 
Brand-new experiments [^5j26| are in progress under the 
current interest. We hope that the second wave of close 
collaboration between theoretical and experimental in- 
vestigations will be brought about. 
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FIG. 1. Temperature dependences of the magnetic sus- 
ceptibility times temperature for the Heiseberg ferrimagnetic 
chains of N = 32. 

FIG. 2. Schematic representations of the M = (S — s)N 
ground state of the decoupled dimers composed of spins 1 
and 1/2 and its ferromagnetic (b) and antiferromagnetic (c) 
excitations. The arrow (the bullet symbol) denotes a spin 1/2 
with its fixed (unfixed) projection value. The solid (broken) 
segment is a singlet (triplet) pair. The circle represents an 
operation of constructing a spin 1 by symmetrizing the two 
spin 1/2's inside. 



FIG. 3. Schematic representations of the M = (S — s)N 
ground states of spin-(S', s) ferrimagnetic chains of N elemen- 
tary cells in the decoupled-dimer limit (a), the AKLT ground 
states of spin-(2s) antiferromagnetic chains of 2N spins (b), 
the M = 2N(S — s) ground states of spin-(S — s) ferromag- 
netic chains of 2iV spins (c). The notation is the same as in 
Fig. 2. The relation (a) « [(b) + (c)]/2 may be expected. 

FIG. 4. Susceptibilities shown in Fig. 1 are compared with 
those of the spin-(S — s) ferromagnetic and spin-(2s) antifer- 
romagnetic Heisenberg chains of N spins. 
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